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Abstract. The freeze out of particles from a layer of finite thickness is dis- 
cussed in a phenomenological kinetic model. The proposed model, based on the 
Modified Boltzman Transport Equation, is Lorentz invariant and can be ap- 
plied equally well for the freeze out layers with space-like and time-like normal 
vectors. It leads to non-equilibrated post freeze out distributions. The depen- 
dence of the resulting distribution on the thickness of the layer is presented 
and discussed for a space-like freeze out scenario. 
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The hydrodynamical modeling of relativistic heavy ion collisions includes three 
main stages: the initial stage, the fluid-dynamical stage and the so-called Freeze 
Out (FO) stage, when the hydrodynamical description breaks down. The FO stage 
is essentially the last part of a collision process and the main source for observables. 

Nowadays, FO is usually simulated in two extreme ways: FO happens on the 
hypersurface with zero thickness, or FO simulated by volume emission model, which 
in principle requires an infinite time and space for a complete FO. For practical 
purposes one has to choose between these two, in some sense opposite, scenarios, 
realized in many different specific models, because it is really very problematic to 
model FO from the first principles. As it was shown recently in Refs. [ 1 during FO 
system has to be described by the Modified Boltzman Transport Equation (Modified 
BTE), which is even more complicated than ordinary BTE. When the characteristic 
length scale, describing the change of the distribution function, becomes smaller 
than mean free path (this always happens at late stages of the FO), then the basic 
assumptions of the BTE get violated, and the expression for the collision integral has 
to be modified to follow the trajectories of colliding particles. One can avoid troubles 
with FO modeling using hydro+cascade two module model [E], since in hadron 
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cascade models gradual FO is realized automatically. In fact, in cascade model one 
follows the trajectories of the colliding particles, therefore what is effectively solved 
is not BTE, but Modified BTE! 

Once the necessity of the Modified BTE and at the same time the difficulty of its 
direct solution were realized, the simplified kinetic FO models become even more 
important for the understanding of the principal features of these phenomenon. 
In the recent papers [ 03 Ej a phcnomcnological kinetic model, which can realize 
complete FO in a FO layer of finite thickness, L, was proposed. This model allows 
us to make a link between the two extreme scenarios, discussed above: L = and 
L — > oo. In fact it is a generalization of a simple kinetic model studied in Refs. [ 
|3J E], which required L — > oo for complete FO. 

According to the above references, the basis of the model, i.e. the invariant 
escape probability within the FO layer of the thickness L, for both time-like and 
space- like normal vectors is given as [0 0] 



L — x^da^J \ p^u 

where p^ 1 is a particle momentum, is the flow velocity and A is some characteristic 
length scale, e.g. mean free path. Here we will concentrate on the space-like case 
only, where the above O function is important. Calculations will be done in the 
reference frame of front where space-like normal vector is da 1 * — (0, 1,0,0). 

The two components of the distribution function, /: the interacting, and 
the frozen out, ft ones, {f — fi + //), evolve according to the following equations 
(for the details of this ID model and its applicability range check Ref. [|3]): 

@fi Pesc p . feqfa) fi ®ff Pesc r /r>\ 

Here the rethermalization of the interacting component is taken into account via 
the relaxation time approximation, where fi approaches the equilibrated Jiittner 
distribution, f eq (x), with a relaxation length, Ao- The system can be solved semi- 
analytically in the immediate rethermalization limit [ 3 . The solution is particularly 
simple for the massless ideal gas without conserved charges. Below we present 
results for such a system aiming for a qualitative illustration of the FO process. 

Figs. IH2l show the results of our calculations. As it was already shown in [0] the 
final post FO particle distributions are non-equilibrated distributions, which deviate 
from thermal ones particularly in the low momentum region. By introducing and 
varying the thickness of the FO layer, we are strongly affecting the evolution of the 
interacting component - see Fig. but for L > 5A the final post FO distribution 
looks very close to our results for an infinitely long FO calculations, see If we 
make the thickness of the layer, L, very small, e.g. L = 2A, then the final post FO 
distribution is practically a thermal distribution. Here FO process does not have 
enough space to distort the thermal shape of the distribution. 

We can propose the following algorithm for a simulation where our idealized 
model would be directly applicable. 1) If the post FO spectra are curved at low 
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Fig. 1. Evolution of the interacting component: upper subplot shows the evolution 
of the temperature of the interacting (equilibrated) component for different FO 
layers, lower subplot - the evolution of its flow velocity. 

pt then FO layer is relatively thick, L > 5A, otherwise we have a narrow FO layer, 
which can be well approximated by FO hypersurface. 2) If the post FO spectrum 
(experimental) is curved, then it doesn't matter how thick was FO layer, we do not 
need to model the FO dynamics in simulations of the collisions. Once we have a 
good parameterization of the post FO spectrum (still asymmetric, non-thermal), it 
is enough to write down the conservation laws and non-decreasing entropy condition 
with this distribution function and probably with some volume scaling factor to 



4 



V.K. Magas et al. 




0.5 1 1.5 

p [GeV] 

Fig. 2. Final post FO distribution for different FO layers as a function of the 

momentum in the FO direction, pt = p x in our case (p y = p z = 0). The initial 

conditions are the same as in Figs. H T(0) = 170 MeV, v(0) — 0.5. The curve for 

L — > oo was calculated in the first version of the model [ 5 . The dotted line shows 

some effective thermal distribution with T = 140 MeV and v = 0.33. 

effectively account for the expansion during FO. Thus, our results may justify the 

use of FO hypersurface in hydrodynamical models for heavy ion collisions, but with 

a proper non-thermal post FO distributions. In large scale simulations for space- like 

FO the Canceling Jiittner distribution [|9j may be a satisfactory approximation. 
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